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By virtue of the generalized Hermann-Feynmam theorem and the method of characteristics we 
rederive energy average formula of photon gas, this is another useful application of the theorem. 
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I. INTRODUCTION 

It is common knowledge in the black-body radiation 
that the expectation value of a Planck oscillator of fre- 
quency vj s is (excluding the zero-point energy term) is 



E, 



i 



(i) 



Historically, Planck derived this result in 1900 first by 
using thermodynamics and then by classical statistics. 
Later in 1924 Bose studied question such as the "proba- 
bility of an energy level E s (= fvuj s ) being occupied by 
n s photons at a time" , he found the mean value of n s is 
eA ihi 3 _ 1 ■ Einstein, after reading Bose's paper, concluded 
that the basic fact to remember during the process of 
distributing the various photons over the various energy 
levels is that the photons are indistinguishable — a fact 
that had been implicitly taken care of in Bose's treatment 
(he essentially treated the light quanta as particles of a 
gas). Planck and Bose-Einstein obtain the same mean 
photon number and mean energy formlas with different 
approaches. We then ask: is there any other method 
with which we can rederive this formula? If yes, then 
this method is worth paying attention in quantum statis- 
tics. The purpose of this work is to re-derive the energy 
average formula of photon gas by using the generalized 
Hermann-Feynman theorem, a theorem regarding to the 
ensemble average of mixed states. 

In quantum mechanics the usual Hermann-Feynman 
theorem [H, Q is for pure states, which states that 



dEj (x) 
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OH (x) 
dx 



li) , U li) = i, 



(2) 



where H is the Hamiltonian depending on the parameter 
X, and Ej and \j) are energy eigenvalues and eigenvec- 
tors of H, respectively. This theorem has been widely 
used in molecular physics, quantum chemistry and quan- 
tum statistics. Having noticed that Eq.([T]) is just for pure 
state expectation value, we need a generalized Hellmann- 
Feynman theorem for mixed state in the sense of ensem- 
ble average in quantum statistics. In Ref. [|[ the general- 
ized Hellmann-Feynman theorem (GHFT) has been set 
up 
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where the subscript e stands for the ensemble average, 
j3 = (fcT) -1 , k is the Boltzmann constant, and (A) = 
Tr (eT 0H A) /Tr (e~ fm ) for arbitrary operator A. When 
H is independent, QM- — 0, Eq.Q can be reformed as 
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Now for the Bose oscillator Hamiltonian H = hzca^a, 
using ([4]) we see 
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that is 
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which can be solved by virtue of the method of charac- 
teristics [H, according to it we have the equation 



dm _ dp _d(H) c 
^~~T~ {H) e 



(7) 



It then follows two integration constants 

lnro = -ln/3 + Cr, ln(H) e -\nm = C 2l (8) 

which tells us that (H) e is of the form 

{H) e = wF{(im) (9) 

the form of F is to be determined. Eq. j9| can be con- 
firmed by noticing 
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dm d In m 
and rewriting ^ as 



, Inn; — x, hxj3 = y, (10) 
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its solution is 



Using the rule of searching for limitation, 



(H) e = e x g(x + y) = mF((3m) 



(12) 



In order to make sure of F's concrete form, we introduce 
P) /V%, so up to a zero-point 



7 -X 



energy term, 
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H = 



2m 



-mm 2 X 2 . 



(13) 



— = lim ,, K = h lira m x -, (20) 

comparing [2^)1 with ([TO)) and considering (fT5|) we can de- 
termine the form of F (f3m) , 



F(J3w) 



1 



(21) 



We see when m — > 0, — * (free particle), Eq. © 
reduces to 



lim , = lim mF (/3m) 



so if lim ro ^o (/^^) is finite, 



lim F (f3zu) — > oo. 



(14) 



(15) 



so follows (H) e = 

Besides, using the GHFT we can also provide: 

1. average kinetic energy = average potential energy. 

To prove this, we notice 



[XPH] =i imm 2 X 2 , 

m 



(22) 



On the other hand, we recall the Bloch equation for the 
density matrix p of a free particle 



dp (P, P) P 2 ( r,. 



its normalized solution 
P(P,P) 
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so the energy average for Hf — 
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Trp = 1. (17) 
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Combining (TJ5]) and fT4]) we see 

(H f ) e = hna (if) e = lim roF (/3tu) = x oo = 



1 

w 



(19) 



2. Though iJ involves the mass, (H) e is m indepen- 
dent. In fact, using (23) and (2) we have 



5m V /e 5/3 



^ \ 5m / 

ski* 



2m 2 



= 0. (25) 



In sum, we have rederived the energy average formula 
of photon gas (Bose distribution) by using the GHFT, 
this is another useful application of the theorem. The 
other application can be seen in Ref . [5j . 
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